Ferromagnetism in one dimension: Critical Temperature 
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Ferromagnetism in one dimension is a novel observation which has been reported in a recent 
work (P. Gambardella et.al., Nature 416, 301 (2002)), and it is thought that anisotropy barriers are 
responsible in that relevant effect. In the present work, transitions between two different magnetic 
ordering phases are obtained as a result of an alternative approach. The critical temperature has 
been estimated by Binder method. Ferromagnetic long range interactions have been included in a 
special Hamiltonian through a power law that decays at large interparticle distance r as r~" for 
Q > 0. We found that if the range of interactions decreases {a oo), the trend of the critical 
temperature disappears, but if the range of interactions increases (a ^ 0), the trend of the critical 
temperature approaches to the mean field approximation. The crossover between two these limiting 
situations is discussed. 

PACS numbers: 02.50.-r; 64.60.-i; 75.10.-b; 75.10.Pq; 75.10.Hk 

Keywords: Thermodynamics, Phase transitions. Magnetic ordering, One dimensional systems. Classical spin 
models 



Recently, much attention has beenpaid to structures 
of lower dimensionality m, |l ll y, |5l Ifl 13 ■ As the space 
dimension of a physical system decreases, magnetic or- 
dering tends to vanish as fluctuations become relatively 
more important. In particular, there are no spontaneous 
magnetization in several one dimensional models, at any 
nonzero temperature; for instance, the isotropic spin- 
S Heisenberg model with finite range exchange interac- 
tion Q and the classical gas model with hard-core and 
finite range interactions!^. However, anomalies such as 
anisotropy properties as microscopic long range interac- 
tions are not taken into account at finite temperature. 

Regardless, ferromagnetism in one dimension has been 
recently reported for monoatomic metal chain of Co con- 
structed on a Pt substrate, with anisotropy barriers 0. 
It was found experimental evidence that the monoatomic 
chains consist of thermally fluctuating segments of ferro- 
magnetically coupled atoms. Chains evolve into a ferro- 
magnetic long range ordered state owing to the presence 
of anisotropy barriers below a threshold temperature . 

Much effort has been devoted to handling finite and 
infinite range interactions in computational systems by 
molecular dynamics and Monte Carlo simulations due to 
the absence of exact and analytical results. However, 
some special situations in one dimension can be studied 
exactly; for instance, 

• mean field theory (this is, if the range of in- 
teractions was infinite) and coupling to first 
neighbors 0, 

• 1/r and logarithmic potentials in a periodic media 
by infinite repetition of a central cellp^. 

• hard core classical gas@, 



isotropic spin-s model 0|, 
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• the thermodynamics of the Casimir force and the 
excess free energy of d-dimensional spherical model 
[13 and 

• a finite size scaling theory is developed when a par- 
ticular family of interactions decays slowly with the 
distance 0, etc. 

In this paper, we present a novel approach to fer- 
romagnetism in one dimension. More explicitly, the 
presence of infinite range microscopic interactions in- 
duces some important modifications to the thermody- 
namic properties of systems, some of them have still been 
not characterized. Several evidences of a ferromagnetic 
state have been suggested due to the existence of long 
range interactions in some physical systems in one di- 
mension II IE 111 111 111 III EE E lla 

The main aim of the present work is to obtain the criti- 
cal temperature between two states of different magnetic 
ordering in a spin-i Ising linear chain, where the range 
of interactions is, at least, comparable to the size of the 
chain. It is well known that there is a state of magnetic 
ordering in the mean field approximation. That approxi- 
mation focusses on a single particle and assumes that the 
most important contribution to the interaction of such 
particle with all particles is determined by the mean field 
due to other particles. The configuration with lowest en- 
ergy is one in which the spins are totally aligned. Before, 
it was explained that no magnetic ordering is observed 
for finite range of interaction {e.g., first neighbors); but, 
which happens between infinite and finite range of inter- 
actions will be illustrated through a generic power law 
decaying as where r is the distance between two 

particles and 0<q;<oo, (q; = and a — s- oo close 
respectively to mean field and to independent spin ap- 
proximations). Hence, the principal question that we try 
to answer in the present work is "How does the critical 
temperature depend on the range of the interaction?" 
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Firstly, the critical temperature for a particular case 
(namely, a = 2) of this kind of systems has been ear- 
lier reported b y va rious authors and by the use of several 
techniques [13, ^ comparison of results is done in a 
previous work, see for instance Ref. 20]. Secondly, the 
critical temperature as a function of other values of a 
(namely, 1 < a < 3/2) was reported in Ref.j^. The 
critical temperature tends to infinite as a tends to 1+ 
(this is, liniQ,^! Tc — > 00) in the thermodynamic limit. In 
general, a strong dependence on the size of the system 
has been observed for such model. Thirdly, the thermo- 
dynamic limit is not reached for a = 1. 

So that, no results have been previously reported for 
a < I. In the present work, we propose a model for 
this kind of system which allows to calculate the critical 
temperature by standard methods in the thermodynamic 
limit for < a < cxd. We carry out a numerical study 
through Monte Carlo procedure and we report direct re- 
sults from our model which considers information on the 
range of interactions. 

The system is described by the following Hamiltonian 



1.5 
1.0 
0.5 
0.0 



0.15 



0.10 



0.05 



0.00 



0.8 



a = 2.0 n = 500 



(a) 



0.2 0.4 0.6 0.8 1.0 1.2 



T/T' 



a = 1.0 


n = 500 


(b) 






0. 
















n = 100~^ 















0.9 



1.0 



1.1 



1.2 



T/T' 



FIG. 1: The susceptibility as a function of the temperature 
T is depicted by (a.) a — 2 and (b) q = 1, the values of n are 
indicated in the figure. 



H = - ^ J{i-j)s^^ 



(1) 



where Si = ±1 Vi and {i — j) is the distance between 
two sites and n is the number of particles in every cell. 



l=L 



1 ^ J tin) 



2 \nl + K\ 



(2) 



where 



(2L-f l)n = N 



is the total number of particles in 2L 
central configuration. 



(3) 

1 repetition of a 



J tin) 
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(4) 



where J is a positive parameter, J^(n) measures the 
strength of the coupling that depends on the size n of 
the system, and J^(n) -> (a - 1) J/2" if a > 1 

Let us consider a computational cell in one dimension 
with size n. Periodic boundary conditions have been ap- 
plied through infinite replications of a central cell. Re- 
cently, the problem related to the periodic boundary con- 
ditions in systems with microscopic long-range interac- 
tions has attracted the attention of several authors (see, 
for instance 20, 21, 23, 24, 25j). Periodic boundary condi- 
tions were applied in a similar manner which was recently 
discussed I2 11. However, this way was already introduced 
by Curilef^. 

Thermodynamics describes the behavior of systems 
with many degrees of freedom after they have reached a 



state of thermal equilibrium. Furthermore, their thermo- 
dynamic state can be specified in terms of a few param- 
eters called state variables. At equilibrium, macroscopic 
observables are linear functions of n (number of parti- 
cles). If the function f{n) is an observable, the thermo- 
dynamic properties impose to observables to be a homo- 
geneous linear function of n, this is f{n) — nf^, where 
fn = fin)/n and /(An) = \f{n) for very large nj^ l2^. 
At this point, it is important to remark that by choosing 
the Hamiltonian of Eq.lQJ (that contains ferromagnetic 
interactions that decay as a l/r" law) we can find two 
facts: 

1. The nice property, about observables as a linear 
homogeneous function of n, is verified for thermo- 
dynamic functions, e.g., the internal energy. 

2. The parameter A can be explicitly written from 
Eq.©. 

As a consequence of the previous properties, a weak de- 
pendence of the size of the system is expected for results 
that come from this model. 

The present computational study was carried out on a 
linear chain with L = 10^ and n between 10^ and 10"^, 
to give effective sizes of the chain of the order from 10^ 
to 10^ particles in according to the Eq.lO, and several 
values of a. Standard techniques are taken into account 
to compute its normalized autocorrelation function C'(t), 
C{t) ^ 1%, due to the large fluctuations near the critical 
point. 

The magnetic susceptibility as a function of tempera- 
ture is estimated from the magnetization fluctuation. 



xiT) 



(5) 



for T>Tc 
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FIG. 2: The specific heat as a function of the temperature T 
is shown. The trend of these functions shows also a disconti- 
nuity. The values of a = 1, 2 and n = 100, 500 are indicated 
in the figure 

where Tc is the critical temperature, < s >= for 
T > Tc- In the FigQ] typical sets of curves for the mag- 
netic susceptibility are depicted as a function of a reduced 
temperature T/T^'^, where T^'^ = J^{n)/kB, where 
fcs is the Boltzmann constant. Range of interactions is 
given by (a) a = 2 and (b) a = 1 , and sizes of the central 
cell are shown for n = 100 and 500, for L = 10^. The 
total number of particles is given by Eq. (|3J . The trend of 
the magnetic susceptibility x(T) suggests a discontinuity 
in Fig^ and a change of the magnetic ordering phase is 
related to such property. Circles are obtained from sim- 
ulations and lines correspond to a linear interpolation of 
the data in the Fig^ It is expected that the trend of the 
magnetic susceptibility becomes a behavior similar to the 
mean field approximation, this is the magnetic suscepti- 
bility has an infinite jump at the critical temperature. 

The specific heat is estimated as a function of the tem- 
perature from the energy fluctuation, which is given by 

C{T)^^{{H\-{H)l). (6) 

In Fig[21 the typical trend of specific heat is depicted. 
Numerical calculations were carried out for several val- 
ues of n and a; namely, n = 100, 500 and a = 1, 2. 
A discontinuity is also observed for the specific heat at 
critical temperature. 

The phase transition can be characterized by several 
ways. A suitable approach to define the critical point in 
finite system is the Binder method. A typical property is 
looked for from the profile of the fourth-order cumulant of 
the magnetization 123. A standard behavior of the phase 
diagram is expected for this kind of system. The Binder 
cumulant of fourth order is defined as 



TABLE I; Comparison of our estimation of the critical tem- 
perature to previous result in the literature ^ij. 



a 


Tc This work 


Tc = {l-a)/2"Kc 


1.1 


0.9689 


0.9797309 


1.2 


0.9625 


0.9438766 


1.3 


0.9101 


0.8951426 


1.4 


0.8813 


0.8367191 


1.5 


0.86 


0.7714285 



Cumulants Un as a function of the temperature, for 
several sizes of system n, are intersected in a common 
point. Such point is the critical temperature which de- 
pends on values of the parameter a. The typical behavior 
of the Binder cumulant Un (or equivalent quantity_o„ = 
— 3C/„ called the renormalized coupling constant[23) is 
rather stable under critical temperature; however, fluc- 
tuations can be important above the critical point. So 
that, this properties permits to distinguish the critical 
point by a very simple graphical criteria. 

In the Figl^l the critical temperature is depicted as 
a function of 1/a for the magnetic ordering transition 
in the thermodynamic limit. Error bars are included to 
represent a 5 percent or less of the deviation for the ob- 
tained values at each point of the critical temperature. 
Additionally, an inset is included in FigO to show the 
same fact for < a < 3. Both pictures are included to 
emphasize the following features about critical tempera- 
ture: 

• It is close to the mean fleld one for a — > 0. 

• It is lesser than the critical temperature given by 
the mean field approximation for a > 1. 

• It falls to zero for the short range interaction regime 
(e.g., nearest neighbor) as it is expected. 

• It is depicted as a function of l/a — > to remark 
that goes to zero while a ^ 00. 

• It is depicted the crossover between two known lim- 
iting cases. 

• It is recovered by the present model the known re- 
sults for particular values of a; namely, a — 1.1, 
1.2 1.3, 1.4 and 1.5'2^ and a = 2^21]. 

We can compare the critical temperature with earlier 
results which a re p resented in the literature. We have 
chosen the Ref.[^ because the authors have made an 
exhaustive search of the critical couplings as a func- 
tion of a and compare their values to others 30). In order 
to compare critical temperatures it is necessary to make 
a simple transformation from Eq.Q for N —^ 00 given 
by Tc = (1 - a)/2°'Kc. The comparison is in the TableU 
In a similar way, for a = 2, our critical temperature is 
0.4221 and we can co mpar e it to a previous value 0.3816 
obtained from the Ref . |20l| . 

In general, according to the Table HI our estimation 
for the critical temperature is greater than the value 
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FIG. 3: The critical temperature as a function of 1/a, the 
inverse of the range of the microscopic interaction, is depicted. 
The trend of Tc goes to zero when 1/a — > 0, namely a — > cx). 
The inset shows the trend of Tc as a function of a. Error bars 
represent a deviation of the 5 percent or less at each point of 
results. 

obtained by other author previously. As it is expected 
that in the limit of short range interactions Tc — » 0, 
it is reasonable to expect that the lower value could be 
more acceptable. Simulations for larger systems had been 
carried out in order to obtain accurate results in one 
dimension _16;j. However, the numerical work is much 
more costly due to the enormous number of particles 
{n = 300000) in systems, in contrast to the method which 
has been implemented here, few particles in a cell and a 
big number of repetitions over all space. If replicas of 
the central cell are on increase, the time of computation 
does not practically suffer changes, because such time 
only depends on the number of particles on the central 
cell. Sometimes, computation facilities are not sufficient 
for filling the demands of resources that make simula- 
tions of many-body systems; then, it is very important 
to resort to an alternative numerical way. 

On the one hand, a simple theoretical argument given 
by L. ReichlQ shows that a finite Ising chain in one 



dimension with a number of ferromagnetically coupled 
spins cannot exhibit a phase transition. If the external 
magnetic field goes to zero the order parameter (s) tends 
also to zero. Hence, no spontaneous nonzero value of the 
order parameter is possible. On the other hand, another 
simple theoretical argument given by the same author 
shows that the mean field theory predicts a phase tran- 
sition at a finite temperature for a lattice in one dimen- 
sion. Both arguments are not opposite between them. 
In the present point of view, the mean field theory (e.g. 
a = 0) approaches to the case of every spin interacts with 
each other without discriminating over the sites and an 
Ising chain is the limit whose interactions are very short 
ranged (e.g. a ~ oo). In this work, we have discussed the 
crossover between both limiting cases, with a system of 
finite ferromagnetically coupled spins that obeys to the 
Hamiltonian given by Eq.lQJ. 

Summarizing, the inhomogeneity which generally char- 
acterizes to systems with interactions with a < 1 is re- 
moved by the scaling J/J^{n) suggested in Eq.lQ, which 
is used to define the Eq.([3l- Thus, if we suppose that 
the scaling represents the number of nearest neighbor 
spins, the size of the system (the total number of spins 
in the chain) is always greater than the such scaling for 
< a < 1. Both values, nearest neighbor and size of 
the system, are coincident for a — 0. In this way, we 
expected that a proper thermodynamic limit can be de- 
fined for a < 1. The scaling is defined and revised by 
several authors (see for instance Ref.Q and references 
therein). With such considerations we can repeat the 
mean field approximation and we hope that the critical 
temperature comes to be exact. In addition, the most 
important problem on the thermodynamic behavior re- 
lated to the systems with long range interactions, it is the 
strong dependence on the size of systems. It is thought 
that no standard thermodynamic equilibrium is reached; 
then, it is crucial to make a suitable choice of the Hamil- 
tonian. However, we have presented a possible way to 
solve such problem. A nice thermodynamic behavior is 
obtained from the Hamiltonian suggested in Eq.J^I to 
Eq.Q). Thermodynamic quantities are not dependent 
on the size of the system for an appropriated choice of 
the Hamiltonian. 

Finally, in the study and characterization of the phase 
transition and critical phenomena for systems with ar- 
bitrary range interactions, advances and suggestions are 
always very important to find appropriated models for 
calculating typical physical quantities. 
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